Comparison of CALIPSO-like, LaRC and MODIS Retrievals of Ice Cloud Properties over 
SIRTA in France and Florida during CRYSTAL-FACE 

M. Chiriaco (1) , H. Chepfer (1) , P. Minnis (2) , M. Haeffelin (1) , S. Platnick (3) , D. Baumgardner <4) , 
Dubuisson (5) , M. McGill <3) , V. Noel (2) , J. Pelon (6) , D. Spangenberg (7) , S. Sun-Mack (8) , G. Wind (9) 

(1) Laboratoire de Meteorologie Dynamique, IPSL, Palaiseau, France 

(2) NASA Langley Research Center, Hampton VA, USA 

(3) NASA Goddard Space Flight Center, Greenbelt, USA 

(4) Universite Autonoma de Mexico, Mexico 

(5) ELICO, Universite du Littoral, France 

(6) Service d’Aeronomie, IPSL, Paris, France 

(7) AS&M, Inc., Hampton, VA, USA 

(8) SAIC, Hampton, VA, USA 

(9) SSAI, Lanham, MD, USA 

Submitted to the Journal of Applied Meteorology 
April 2005 

Corresponding author: 

Marjolaine Chiriaco 

33 (0)1 69 33 31 94 

chiriaco@lmd.polvtechnique.fr 

Laboratoire de Meteorologie Dynamique / IPSL 

Ecole Polytechnique 


91128 Palaiseau cedex, France 


Abstract 


This study compares cirrus particle effective radius retrieved by a CALIPSO-like method with two 
similar methods using MODIS, MODIS Airborne Simulator (MAS), and GOES imagery. The 
CALIPSO-like method uses lidar measurements coupled with the split-window technique that uses 
the infrared spectral infonnation contained at the 8.65-pm, I 1 . 1 5 -pm and 1 2.05-pm bands to infer 
the microphysical properties of cirrus clouds. The two other methods, using passive remote sensing 
at visible and infrared wavelengths, are the operational MODIS cloud products (referred to by its 
archival product identifier MOD06 for MODIS Terra) and MODIS retrievals performed by the 
CERES team at LaRC (Langley Research Center) in support of CERES algorithms; the two 
algorithms will be referred to as MOD06- and LaRC-method, respectively. 

The three techniques are compared at two different latitudes : (i) the mid-latitude ice clouds study 
uses 18 days of observations at the Palaiseau ground-based site in France (SIRTA: Site Instrumental 
de Recherche par Teledetection Atmospherique) including a ground-based 532-nm lidar and the 
Moderate resolution Imaging Spectrometer (MODIS) overpasses on the Terra Platform, (ii) the 
tropical ice clouds study uses 14 different flight legs of observations collected in Florida, during the 
intensive field experiment CRYSTAL-FACE (Cirrus Regional Study of Tropical Anvils and cirrus 
Layers - Florida Area Cirrus Experiment), including the airborne Cloud Physics Lidar (CPL) and the 
MAS. 

The comparison of the three methods gives consistent results for the particle effective radius and the 
optical thickness, but discrepancies in cloud detection and altitudes. The study confirms the value of 
an active remote-sensing method (CALIPSO-like) for the study of sub-visible ice clouds, in both 
mid-latitudes and tropics. Nevertheless, this method is not reliable in optically very thick tropical ice 
clouds. 
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1 . 


Introduction 


Ice clouds play a major role in the radiative energy budget of the earth-atmosphere system 
(Liou 1986). Their impact is governed primarily by their opposing albedo and greenhouse effects that 
are poorly known for thin ice clouds. Both the macro- and microphysical properties of ice clouds 
detennine these competing effects. The global coverage, altitude, temperature, vertical structure and 
spatial inhomogeneities of these clouds must be accurately assessed to quantify their impact on 
weather and climate. Additionally, ice water content and its spatial distribution are critical to the 
global radiative effect of cirrus clouds. A large uncertainty in evaluating the radiative impact of ice 
clouds arises from our limited knowledge of the natural variability of their microphysical properties, 
such as the size, shape and orientation of ice crystals, which detennine their optical characteristics. 
For example, the effective radius of ice crystals composing cirrus clouds is an important 
microphysical parameter, as it largely affects cloud albedo, which increases with decreasing effective 
radius for a fixed ice water path. Reducing the uncertainties in both their macro- and microphysical 
properties is essential for accurately defining the role of ice clouds in climate. 

Satellite data are best suited for estimating ice cloud properties because of their global 
coverage and long-term monitoring potential. Many remote sensing techniques have been devised to 
take advantage of the spectral radiances measured from various satellites to retrieve ice cloud 
properties. The remote sensing methods have been complemented on a local scale by various mid- 
latitude field experiments, such as the First ISCCP Regional Experiment (FIRE; Randall et ah, 
1996), the European Cloud and Radiation Experiment (Raschke et ah, 1998), the International Cirrus 
Experiment (Raschke et ah, 1990), and the SUbsonic aircraft: Contrail and Cloud Effects Special 
Study (SUCCESS; Toon and Miake-Lye, 1998) as well as a few tropical campaigns such as the 
Central Equatorial Pacific Experiment (CEPEX; McFarquhar and Heymsfield, 1996) and the recent 
Cirrus Regional Study of Tropical Anvils and cirrus Layers - Florida Area Cirrus Experiment 
(CRYSTAL-FACE; see Jensen et ah, 2004) have collected a wealth of in-situ data that show 
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significant variability of cirrus cloud particle sizes, either from one cloud to another, or within a 
single cloud (Heymsfield 1975, Heymsfield and Platt 1984, Heymsfield 1993, Krupp 1991). 
Improved knowledge of realistic effective radii and shapes of ice crystals on the basis of in-situ 
measurement has led to an improvement in the computation of the scattering properties of ice 
crystals (Wendling et ah, 1979, Takano and Liou 1989, Macke et ah, 1996, Yang and Liou 1998). In 
parallel, various remote sensing retrieval techniques have been developed to infer the ice crystal 
effective radii from ground-based observations using active sensors (Platt et ah, 1989, Intrieri et ah, 
1993, Intrieri et ah, 1995, Mace et ah, 1998, Matrosov 1999) and from satellite observations using 
passive sensors (Minnis et ah, 1998, King et ah, 2003, Platnick et al., 2003). In addition to improving 
the theoretical treatment of cloud optical properties, the in situ data are useful for understanding and 
evaluating the ground-based and, to some extent, the satellite retrievals. The ground-based retrievals 
can also be used for verifying the satellite passive remote sensing (e.g., Mace et al., 2004). It is also 
necessary, however, to determine the level of agreement between different remote sensing methods 
as another aspect of the uncertainties in the retrieved cloud properties. 

A new retrieval algorithm recently developed for the upcoming Cloud-Aerosol Lidar and 
Infrared Pathfinder Satellite Observations (CALIPSO) mission (Winker et al., 2003) combines both 
passive and active sensors to retrieve ice particle effective radius. This hybrid technique (Chiriaco et 
al., 2004) will be implemented to create an operationally archived CALIPSO product. It will use 
collocated data from the three-channel (8.7, 10.5, and 12 pm) CALIPSO Infrared Imager (HR) and 
the 532-nm backscatter lidar operating with polarization capability on CALIPSO. Given the expected 
CALIPSO launch in 2005, it is important to ensure that the technique operates as expected. The 
intent of this paper is to perform an initial validation of the ice particle effective radius retrieval 
algorithm developed for this instrument by Chiriaco et al. (2004) and to compare a few other relevant 
parameters such as optical thickness and cloud height. This new algorithm is applied to datasets of 
mid-latitude ice clouds observed at the Site Instrumental de Recherche par Teledetection 
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Atmospherique (SIRTA, Haeffelin et al. 2005) in France and of tropical ice clouds observed during 
the 2002 CRYSTAL-FACE field campaign. The results of these retrievals are compared to two other 
methods that use passive remote-sensing measurements from the MODerate resolution Imaging 
Spectroradiometer (MODIS) imager: the algorithm developed for the Clouds and Earth’s Radiant 
Energy System (CERES) measurements by Minnis et al (1995, 1998), and the algorithm developed 
as part of the operational MODIS cloud product (Platnick et al., 2003). The former method will be 
referred to as the LaRC algorithm, while the latter is referred to as MOD06 which is the product 
identifier for the MODIS Terra cloud product in the NASA Goddard Earth Sciences Distributed 
Active Archive Center (DAAC). LaRC retrievals for a CERES footprint are archived at the NASA 
Langley DAAC. Differences between MOD06 and a similar algorithm applied to the MODIS 
Airborne Simulator (MAS) data obtained during CRYSTAL-FACE will be discussed in section 4. 
For CRYSTAL-FACE analysis, LARC retrievals are applied to GOES imagery. In situ data taken 
during CRYSTAL-FACE are also used to evaluate both the passive and CALIPSO-like algorithms. 


2. Three methods to retrieve ice cloud properties 
a. CALIPSO-like method 

The CALIPSO-like method retrieves ice cloud properties, in particular, the optical thickness 
and the effective radius of the ice crystals given a set of measurements similar to those that will be 
taken by CALIPSO: 532-nm lidar backscattering and brightness temperatures at 8.7, 10.5, and 12 
pm. The effective radius r e initially corresponds to the radius of a sphere with a volume equivalent to 
the ice crystal one, but it is corrected in order to be defined as: 

3 V 

(1) 

where V is the particle volume and A is its projected area. 

The method (Chiriaco et al., 2004) is based on the Brightness Temperature Differences (BTD) 
between pairs of infrared channels measured from space. The BTDs are compared to their 
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counterparts computed for various ice crystal models based on effective radius and aspect ratio 
derived by Yang et al. (2001) with optical thicknesses ranging from 0 to 50. The brightness 
temperatures are computed using a radiative transfer code including both scattering and absorption 
(Dubuisson et al., 2005). The cloud base and top altitudes are detennined with the 532-nm lidar. 
Those altitudes are associated with a collocated radiosonde temperature profile to determine the 
cloud top and base temperatures. Then, for each cloud case studied, the surface and clear-sky 
atmospheric properties, as well as the observed cloud temperature at the altitude where the lidar 
backscattered signal is maximal (Chiriaco et al. in revision), are prescribed in the radiative transfer 
code. This method is a split window technique, improved by the lidar information. This improvement 
strongly decreases the uncertainty on effective radius retrieval for semi-transparent clouds (Chiriaco 
et al. 2004). 

For a cloudy atmosphere, various microphysical parameters of particles are used to specify 
the cloud optical properties in the radiative transfer computations: spherical models and hexagonal 
models with aspect ratios varying between 0.15 and 9.26, for effective radii varying between 0.5 and 
200 pm (Chiriaco et al., 2004). The comparison between the computed and measured BTDs is used 
to select pairs of r e and particle aspect ratio Q for several microphysical models that best match the 
observations. Another constraint uses Q derived from the lidar depolarization with the method of 
Noel et al. (2002). Four Q classes are defined to roughly segregate particles from the most to the 
least spherical: class I for Q < 0.05, class II for 0.05 < Q < 0.7, class III for 0.7 < Q < 1.1, class IV 
for Q > 1.1. This constraint reduces the number of solutions, leading to an optimized selection of a 
microphysical model as defined by effective radius and particle aspect ratio. 

b. LaRC method 

The LaRC method was developed for application to operational meteorological satellite 
imagers for the Atmospheric Radiation Measurement program (ARM, Ackerman and Stokes, 2003) 
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and for the research satellite imagers used by the CERES project. This method detects clouds and 
retrieves cloud properties using radiances measured by at 0.65, 3.75, 10.8, and 12.05 pm; effective 
particle radius information comes from the 3.75 pm band, but its definition is different from the r e 
definition in (1). In order to compare effective radii from the same definition, we use this correction 
formula: 

r e = 0.0025 x rl + 0.605 x r 0 (2) 

where ro is the initial radius from LaRC technique. 

The LaRC scene classification technique consists primarily of cascading threshold tests 
(Trepte et al., 2002). To define a pixel as cloudy, at least one of the spectral radiances must differ 
significantly from the corresponding expected clear-sky radiances. For a pixel classified as cloudy, 
the main algorithm used to retrieve cloud properties is the Visible Infrared Solar-Infrared Split- 
window Technique (VISST), which is a 4-channel upgrade of the 3-channel method of Minnis et al. 
(1995). Given the clear-sky radiances and surface properties, the VISST computes the spectral 
radiances expected for both liquid-droplet and ice-water clouds for a range of optical thickness x = 
0.25 to 128 for a particular effective cloud temperature T c using a set of reflectance and emittance 
parameterizations (Minnis et al., 1998). The liquid particle effective radii r e ,iiq U id for the model clouds 
range from 2 to 32 pm and the effective radii r e .i cc for the hexagonal ice column model clouds vary 
from 3 to 67.5 pm. The thennodynamical phase is determined by a combination of tests that 
incorporate the final cloud temperature, the initially derived cloud altitude, and the 10.8 - 12.0 pm 
brightness temperature difference. When the spectral data are available, the 1.6-pm/0.65-pm 
reflectance ratio is also used in the phase decision. The effective cloud height is estimated from T c 
using vertical profiles of temperature from numerical weather analyses. Cloud-top and base heights 
are then estimated from empirical data (e.g., Chakrapani et al., 2002). The LaRC technique has been 
applied to a variety of datasets including the Geostationary Operational Environmental Satellite 
(GOES) series (e.g., Minnis et al., 2004), the National Oceanic and Atmospheric Administration / 
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Advanced Very High Resolution Radiometers (NOAA/AVHRR, Young et al., 1998), the Visible 
Infrared Scanner on the Tropical Rainfall Measuring Mission (TRMM) satellite, and the MODIS on 
the Terra and Aqua satellites (Minnis et ah, 2003). As shown in Table 1, the same version of LaRC 
method is used for both mid-latitude and tropical ice clouds. 


c. MODIS operational method: MOD06 

The MOD06 uses several channels of the instrument ranging between 0.65 and 13.3 pm 
(Platnick et al., 2003; King et al., 2003). The cloud mask uses a series of threshold tests to detect the 
presence of clouds with 4 confidence levels, with up to 20 bands from visible through the infrared 
(Ackerman et al., 1998). For mid- to high-level clouds, the cloud top pressure is retrieved using a 
CO 2 slicing technique (Menzel et al., 1983) that uses ratios of differences in radiances between 
cloudy and clear-sky regions at two nearby wavelengths. For low-level clouds, the infrared window 
11 -pm band temperature is used to determine a cloud top temperature, assuming the cloud is 
optically thick, and a cloud top pressure is assigned by comparing the measured brightness 
temperature to the National Center for Environmental Prediction (NCEP) Global Data Assimilation 
System (GDAS) temperature profile (Derber et al., 1991). For the thennodynamic phase retrieval, the 
method uses the 8.52- and 11-pm bands to exploit the spectral emission differences between water 
droplets and ice crystals (Baum et al., 2000) as well as ratios of the 1.6 and 2.1 pm bands relative to 
visible and/or near-infrared bands (King et al., 2004). Cloud optical thickness and r e are estimated 
based on library calculations of plane-parallel homogeneous clouds overlying a black surface in the 
absence of atmosphere. Separate libraries exist for liquid water and ice clouds, the latter consisting of 
12 size distributions composed of 4 habits (aggregates, bullet rosettes, hollow columns, and plates) 
with the fraction of each habit in individual size distribution bins being a function of particle 
effective radius that is defined as follow: 


3 <V> 

4 < A > 


( 3 ) 
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This definition is equivalent to (1) but for crystals with size distribution: <V> is the mean particle 
volume for the ice crystal size distribution and <A> is the mean projected area. 

Scattering calculations are made using the techniques of Yang and Liou (1996). The effective radii 
for these 12 distributions vary from about 6-60 pm (King et ah, 2004). The optical thickness and 
effective radius retrievals are derived from MODIS water-absorbing bands (1.6, 2.1, 3.7 pm) in 
conjunction with one of the non-absorbing bands (0.65, 0.86, 1.2 pm). 

Table 1 shows the differences between the products available for the SIRTA cases and the 
CRYSTAL-FACE cases. The MOD06 retrievals over SIRTA use MODIS/Terra data and are 
archived at the NASA Goddard Distributed Active Archive Center (DAAC). The CRYSTAL-FACE 
cases are processed independently using radiances from the MAS that flew on the NASA ER-2 
aircraft. The MAS retrievals are based on the retrieval method of Platnick et al. (2001) and use ice 
libraries similar to those described above for MOD06, but with two distinctions. First, new scattering 
calculations for the 12 distributions were made using MAS-specific spectral band-passes. Second, a 
quadratic polynomial was fitted to all scattering parameters (single scattering, albedo, asymmetry 
parameter, and extinction efficiency) as a function of the 12 size distributions in order to smooth out 
non-monotonic behavior with particle effective radius. Then reflectance libraries were calculated at 
equal intervals in interpolated r e space (5, 10... 55, 60 pm). The MAS products include a map of 
apparent multi-layer clouds that provides information on the confidence of the retrievals. Being a 
spectrometer, the phase algorithm for MAS includes additional tests which are derived from the 
location of the MAS band corresponding to the peak reflectance in the 1.6 and 2.1 pm spectral 
regions (Pilewskie and Twomey, 1987). Finally, the MAS cloud top properties retrievals are not 
derived from a CO 2 slicing technique, but from an algorithm that uses path absorption in the 0.94-pm 
water vapor band; NCEP profiles are used to convert vapor path to cloud-top pressure and 
temperature. 
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3. Mid-latitude ice clouds: SIRTA cases 


a. Observations: instruments and dates 

The SIRTA observatory sits on a 10-km plateau about 160 m above sea level. The plateau is a 
semi-urban environment divided equally in agricultural fields, wooded areas, and housing and 
industrial developments. Eighteen mid-latitude cloud cases observed at the SIRTA ground-based site 
are investigated using observations collected with the instruments listed below. 

1) MODIS. Depending on the orbit, MODIS on Terra observes Palaiseau once or twice per 
day. Although higher resolution data are available for selected channels, the nominal MODIS 
resolution is 1 km for all 36 channels. Those data are used in each of the three methods and 
the pixel areas used for the application are described in Table 2. For all retrievals except 
MOD06 cloud top pressure and temperature, the sample of pixels is a strip large of 3 pixels 
wide (i.e. 3 km) and 1 hour long centered on the SIRTA overpass time (using wind 
information from radiosonde at the cloud altitude from lidar); for this sample, data are 
corrected from parallax using the cloud altitude from lidar and the wind speed from the 
radiosonde. For the MOD06 pressure and temperature retrievals, the sample used is a box of 
5x5 pixels (i.e. 5x5 km 2 ) over SIRTA, as this product is only available at this resolution. 

2) Ground-based 532-nrn lidar. The SIRTA 532-nm lidar is similar to the one that will be on- 
board the CALIPSO platform. It operates 4 days per week from 0800-2000 hours local time, 
with a nominal temporal resolution of 10 seconds (an average of 200 shots) and a vertical 
resolution of 15 m. It is a zenith viewing lidar that measures both the backscattered signal and 
linear depolarization ratio. The lidar observations are treated as follows: (i) the parallel- 
polarized signal at 532-nm is normalized to the molecular signal deduced from the radiosonde 
profiles (described below), (ii) the cloud base and top altitudes are derived from lidar profiles 
averaged over 1 minute in the parallel channel using a threshold method, (iii) the linear 
depolarization ratio is the ratio of the perpendicular to the parallel channel and it is 
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normalized to 2.8% (Young, 1980) in a cloud-free area, i.e. where there are only atmospheric 
molecules. 

3) MeteoFrance radiosondes. The ground-based cloud-top and -base temperatures are derived 
from the coupling of the cloud altitude (derived from the lidar) with the temperature profile 
obtained by radiosondes launched twice per day (0000 UTC and 1200 UTC) from the 
Trappes station as part of the Meteo-France operational network, located 15 km away from 
the SIRTA site. 

The studied cases were selected using lidar information. Ice clouds were identified using 
the depolarization ratio based on the behavior of non-spherical particles (i.e., ice crystals) that 
strongly depolarize the signal. Furthermore, the boundaries of thin ice clouds are detennined from 
the lidar backscattered signal when both cloud bottom and top are evident in the returns. Only those 
clouds that are relatively constant (+/- 30 minutes around the satellite overpass) in time and in space 
are selected in order to facilitate good collocation between the lidar and the satellite measurements. 
A total of 18 different days of observations have been analyzed. Table 3 summarizes the dates and 
time periods of observations for the different instruments as well as the cloud altitudes. 

b. Comparison of the three methods 
A) Illustrative case: 8 October 2002 

Figure la shows the evolution of the 532-nm lidar depolarization ratio with time and altitude. 
This image indicates a cloud located between 6.5 and 11 km with a relatively constant altitude, 
making it a good case study. The depolarization ratios range between 20% and 60% within the cloud, 
indicating that it is an ice cloud. Figure lb shows a pseudo-color image of the area based on the 0.6, 
1.6, and 11-pm MODIS channels. Pixels from this image were collocated with the lidar 
measurements corresponding to time when Terra passed over the SIRTA site (center of the red 
square) at 1105 UTC. The length of the strip of pixels used for the retrieval is calculated as follow: 
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each pixel is at a distance d from the SIRTA pixel; the horizontal wind (U, V) at the cloud altitude is 
known from the radiosonde; and the time when the pixel overpasses the SIRTA is calculated from d 
and (U, V). The pixel strip selected for the retrieval is 1 km wide and corresponds to +/- 30 min 
centered on the SIRTA site. This strip is represented with red lines in Figure la. For this case, U = 
5.79 m.s' 1 and V = 2.43 m.s' 1 at the altitude of cloud center, so the strip length is approximately 22 
km. Figure lc shows the cloud phase retrieved from the LaRC method. In the red square centered on 
the SIRTA site, most of pixels are ice clouds, with a few pixels of water cloud or clear sky near the 
cloud edges. 

Table 4 shows the results of the retrievals using the three methods. The MOD06 cloud top 
altitude is derived from the MOD06 cloud top pressure and the radiosonde. The three methods give 
consistent altitudes and temperatures, even if MOD06 and LaRC techniques show a cloud top above 
the lidar one (respectively 1.3 and 2.3 km higher). The mean effective cloud altitude given by the 
LaRC method is 0.8 km above the lidar one (CALIPSO-like). The effective temperature refers to the 
radiative equivalent temperature of the cloud and is typically below cloud center for optically thin ice 
clouds. The visible optical thicknesses obtained with the three methods are consistent, ranging 
between 1.3 and 3. For such a semi-transparent cloud, the CALIPSO-like method uses the lidar- 
based optical thickness, which is corrected for multiple scattering following Platt (1973) using r|=0.5 
(where T] is the ratio of apparent visible extinction coefficient to true visible extinction coefficient). 
The particle effective radii are also consistent varying between, but the CALIPSO-like method gives 
a smaller result. The MOD06 and CALIPSO-like methods give a larger variability (respectively 23 
and 16 pm ) than the LaRC technique (9 pm). 

B) All SIRTA cases 

All the 18 SIRTA cases are studied as the 8 October 2002 case. 

Before any interpretation, it is necessary to note that the MOD06 altitudes and temperatures 
given in Table 3 do not use the same sample of pixels than all the other products. One consequence is 
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that for 5 cases, MOD06 technique does not detect any ice cloud using the 5x5 pixels sample (the 
one used for altitude and temperature retrievals) whereas there is an ice cloud using the pixels strip 
sample (the one used for the other retrievals: optical depth, effective radius, liquid/ice water path). 
Figure 2 is the MOD06 optical thickness retrieved from the 5x5 pixels sample, against the one 
retrieved from the pixels strip sample, when both lead to ice cloud. For 3 cases out of 6, the result is 
different and illustrates the variability that must be taken into account when we compare cloud 
temperature and altitude from the different methods. 

For two cases (2 April and 14 October 2002) in Table 3, the LaRC and MOD06 effective 
cloud altitudes are respectively similar and greater than the cloud top detected by the lidar, with 
consistent clout top temperatures. For those cases, the optical thicknesses determined with the LaRC 
and MOD06 methods are greater than the lidar values. In particular, for the 14 October case, the 
LaRC and MOD06 methods yield Xi aRC/Monoe = 1 1 and 1 1.5 respectively, compared to only 1-3 from 
the lidar. For those large values of optical thickness, the lidar cannot detect the cloud top. Hence, the 
CALIPSO-like particle size retrieval is not reliable for those 2 cases. This limitation of CALIPSO- 
like method to cases where x < 3 is important for ground base applications but will not be as crucial 
for future spatial observations, as both lidar and infrared radiometer will be down-looking allowing 
for the studying of upper cloud layer properties. 

In all cases, the CALIPSO-like method detects ice clouds only. For six cases in 2003 (19 and 
24 February, 27 March, 15 September, 6 November and 9 December), the cloud is too thin for the 
MOD06 method to detect clouds along the wind strips. The same thing is observed for 4 cases using 
LaRC method (19 February, 27 March, 15 September, and 6 November). Those days were selected 
specifically because the lidar detects subvisible ice clouds with optical thicknesses smaller than 0.2. 
The 10 remaining cases are used to compare cloud altitude, optical thickness and effective radius. 
The clouds are identified as mixed water and ice by LaRC method in nine instances, and by MOD06 
in three instances, whereas the clouds are only ice from the CALIPSO-like method. Examination of 
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the imagery reveals that low level clouds are present in the vicinity for the 1 st and 8 October and 
could easily have been included in the pixel strips because of the parallax and advection corrections, 
which would not be properly applied for low clouds. No low clouds were evident for the 9 December 
case, which mostly consisted of contrails. For the 6 November 2002 case, the LaRC method retrieves 
both ice and water clouds. The imagery (not shown) reveals the presence of an extensive, broken low 
cloud deck. The “comments” column of Table 3 indicates when the LaRC and MOD06 methods find 
several ice or water pixels, meaning that the pixels composing the strip are classified as ice and 
liquid phases (a single pixel being associated to one phase only). 

Comparisons for the ten scenes where the lidar optical thickness exceeds 0.3 are shown in 
Figs. 2-5. Figure 3a shows the scatter plot of cloud temperature retrieved with CALIPSO-like 
method (T C AL-iike), as a function of the LaRC retrieved temperature method (T LaRC ). In the graph, the 
variability is the difference between the cloud top and bottom, and the LaRC cloud effective 
temperature is the cross on errorbars. Figure 3d is the probability density function of the difference 
between CALIPSO-like and LaRC cloud temperatures. Those two figures show that T CA L-iike is 
always less than T LaRC , even if they are very close for cases when the CALIPSO-like optical 
thickness (XcAL-iike) ranges from 0.8 to 1. Figure 3d shows that for 75% of cases, the difference 
between TcAL-iike and T LaR c is larger than 20 K and for 25%, the difference is between 0 and 20 K. 
The greatest difference, 35 K, occurred 9 September 2003 associated with a sub-visible cloud with 
XcAL-iike ranging between 0.09 and 0.1 (Table 3). Figure 3b is the comparison between the CALIPSO- 
like temperatures and the MOD06 ones when MOD06 using the 5x5 pixels sample finds an ice 
cloud. As the MOD06 are only available concerning the cloud top, their variability are a standard 
deviation calculated in the sample. We must compare them with the lowest point of CALIPSO-like 
errorbars: it shows that the agreement is very good, even if TcAL-iike is most of the time (70%) greater 
than T M odo6 as shown by Figure 3e. The comparison between T M odo6 and T LaR c (Figure 3c and f) is 
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consistent with the other ones: T LaRC is always wanner than T M odo6, showing than the C02 slicing 
technique used for the MOD06 method is more perfonning. 

The results for altitudes (Figure 4) are consistent with the ones obtained for temperatures: for 
most of the cases, z C AL-iike is greater than z LaRC , except for one case, in which T C AL-iike = T LaRC (5 
March 2002). The greatest difference between z C AL-iike and z LaRC is 6.5 km and corresponds to 17 
November 2003, when XcAL-iike varied from 0.05 to 1. As for the temperature, the agreement is better 
between CALIPSO-like and MOD06 methods. 

Figure 5 a shows the variations of T C AL-iike - T LaR c as a function of XcAL-iike- T M odo6 and z M odo6 
are not considered here because of the too few cases where MOD06 detects an ice cloud using the 
5x5 pixels sample. The difference between the temperatures and altitudes (Figure 5b) generally 
decreases as XcAL-iike increases. The vertical variability is the difference between the cloud top and 
bottom from CALIPSO-like method added to standard deviation from LaRC method. 

Figure 6 shows the comparison between the optical thicknesses without the two cases with 
clouds too thick for the CALIPSO-like retrieval. All the other cases are considered in these statistics, 
even the ones where LaRC and MOD06 methods do not detect any cloud, i.e. when x LaR c/MO06 = 0. 
Figure 6a, b and c show that there is good agreement between the 3 methods for the cases that are not 
too thick for the CALIPSO-like method and not too thin for LaRC and MOD06 methods, typically 
meaning x = 0.5-3. Nevertheless, for 4 thin cases, XcAL-iike is smaller than x LaR c and x M odo6- This is 
certainly due to a collocation problem: XcAL-iike is calculated using lidar so only for the point above 
the SIRTA; x LaR c and x M odo6 are calculated using the strip of pixels that could include a few pixels of 
a thick cloud, leading to larger values of optical thicknesses. Figure 6d shows that the difference 
between XcAL-iike and Xl 3 rc is less than 0.2 for 50% of the cases, a result similar to the difference 
between XcAL-iike and x M odo6 (Figure 6e). Nevertheless, x C AL-iike is smaller than x M odo6 in 45% of the 
cases. There is also a good agreement between Xl 3 rc and Xmodo6 in 45% of the cases (difference 
smaller than 0.2). Overall, Xmodo6 tends to be larger than Xl 3 rc- 
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Figure 7 summarizes the results for the ice particle effective radius. Figure 7a, b and c show 
that while there is generally good agreement between the 3 methods, the variability of CALIPSO-like 
effective radius (r e ,cAL-iike) is often larger than that for either LaRC (r c l aRC ) or MOD06 (r e ,MOD06)- The 
variability difference is probably a result of the differences in the microphysical models used for the 
retrievals. The CALIPSO-like and LaRC methods have the best agreement: in 60% of the cases, the 
difference between mean values of r ejC AL-iike and r ej LaRC is less than 3 pm (2 main peaks). In 35% of 
the cases, r e>C AL-iike is at least 5 pm larger than r c laRC . This difference can be very important, in 
particular for 5 March 2002 case where one can only know that r e cAL-iike is greater than 14 pm. The 
comparison between r ejC AL-iike and r ejM OD06 (Figure 7b and e) is similar, but the differences are slightly 
greater: for 45% of the cases, r ejC AL-iike is at least 5 pm smaller than r ejM OD06- For 85% of cases, 
r e ,MOD06 larger than r e!LaRC and this difference is at least 5 pm for 75% of cases (Figure 7c). 

4. Tropical ice clouds during CRYSTAL-FACE field experiment 
a. Observations: instruments and dates 

Ice cloud observations collected during CRYSTAL-FACE in July 2002 are used to evaluate 
the CALIPSO-like method in the Tropics. 

1) Airborne 532-nrn lidar. The CPL (Cloud Physics Lidar) is a Nd: YV04 laser. It was nadir- 
viewing onboard the high altitude NASA ER-2 aircraft, and measures both the backscattered 
signal at 532 nm and the linear depolarization ratio at 1064 mn. The temporal resolution is 1 
Hz, and its vertical resolution is 30 m. For the current study, cloud-base and -top altitudes and 
the number of cloud layers are from the CPL archived files (http://cpl.gsfc.nasa.gov/). The 
CPL cloud temperatures were found by comparing the cloud base and height data with 
vertical profiles from the ER-2 dropsondes. The crystal shape information is based on 
analysis of the lidar depolarization data following the methodology of Noel et al. (2002). 
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2) Satellite-borne imager. The 4-km GOES-8 imager provided observations every 10-15 
minutes enabling relatively close matches with the aircraft data. The GOES-8 data were 
analyzed with the LaRC method (Table 2) using hourly temperature profiles from the Rapid 
Update Cycle analyses. The retrieved pixel-level cloud properties were averaged for groups 
of the 4 closest pixels to a point on ER-2 flight track. 

3) An airborne spectrometer. The MAS instrument flown on the NASA ER-2 aircraft has 
similar capability as the MODIS instrument; although with much higher spatial resolution (50 
m) and more complete spectral coverage provided by grating spectrometers (50 spectral 
bands). The MAS data are used for both the MAS retrieval and CALIPSO-like methods using 
a strip of pixels that is 3 pixels (150 m) in the across-track direction, along the length of the 
selected flight leg. The cloud properties are retrieved for each pixel and then averaged 
together over the flight leg. The MAS retrieval has similarities to the MOD06 retrieval as 
previously described, especially in its use of the 2.1 pm band for the primary particle size 
retrieval. 

4) ER-2 Dropsondes. These dropsondes are an automated version of the NCAR Airborne 
Vertical Atmospheric Profiling System (AAV APS) produced by Vaisala and launched from 
the ER-2. Details pertaining to the AAV APS characteristics and performance are described in 
Hock and Franklin (1999). The AAV APS provides measurements of pressure, temperature, 
humidity, and horizontal wind. The sampling rates are 0.5 and 1 Hz for thermodynamic and 
wind data, respectively. From an ejection altitude of 24 km, the dropsonde requires about 22 
minutes to descend to the ocean surface. Between 5 and 6 dropsondes were released during 
each ER-2 flight. Dropsondes are not available for 29 July flight; hence for this day, data 
taken by the Balloon Borne Sounding System (BBSS) launched from the Key West ground 
station (25° N, 81° W) were used. 
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5) In-situ measurements . Hydrometeor size spectra built from combined datasets taken by the 
CAPS (Cloud, Aerosols, and Precipitation Spectrometer), CPI (Cloud Particle Imager), and 
SPP-100 optical spectrometer on the WB-57 aircraft are used to evaluate the remote sensing 
techniques. The CAPS measures ice particle diameter over two ranges, 0.5 - 44 pm and 75 - 
1600 pm (Baumgardner et ah, 2001); the SPP-100 measures between 4 and 60 pm, while the 
CPI covers 10 - 300 pm. After re -binning each measurement, an average concentration is 
computed for each overlapping effective size bin. The particle volume and cross sectional 
area were computed from the size distributions with the assumption that the particles were ice 
spheres with density 0.9. Both parameters will be overestimated since recent studies 
(Baumgardner et ah, 2005) show that the cirrus particles in CRYSTAL-FACE had shapes that 
are representative of ensembles of bullet rosettes, plates and hollow columns. All of these 
habits have equivalent volumes less than a sphere. The uncertainty in volume is 
approximately ±50%, based on comparisons of different habits with equivalent spheres. We 
estimate the uncertainty in area as approximately ±30% since the cross sectional areas of 
rosettes, columns and plates are much more similar to those of spheres with the same 
equivalent diameter. The uncertainties in volume and area are correlated; hence, error 
propagated into the calculation of effective radius is approximately ±40%. 

Collocated ER-2 and WB-57 aircraft data were determined using navigational recorder data 
provided by K. Drdla (also available at http://www.espoarchive.nasa.gov/cgi-bin/dl_start). 
The ER-2 and WB-57 data are considered to be coincident when the planes pass the same 
location within + 5 minutes. 

The cases were selected using almost the same criteria as in mid-latitudes (Sect. 3.b). The 
main difference is that for CRYSTAL-FACE, each case corresponds to a flight leg instead of a day, 
and several flight legs are studied in a given day. Using the CPL information, a flight leg is selected 
when the cloud has relatively constant (+/- 150 m) bottom and top altitudes so that a single IR 
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radiative transfer computation can be used for the entire flight leg to apply the CALIPSO-like 
method. The lidar depolarization is used to select ice cloud in flight legs. For those cases, the 
CALIPSO-like method is not limited to cases where x < 3 because both the imager (MAS or GOES) 
and lidar (CPL) are down-viewing. Nevertheless, other limitations of this method will be 
demonstrated. A total of 14 different flight legs were analyzed. Particle effective radii and optical 
thickness were retrieved for each pixel and then averaged along the flight leg. Table 5 summarizes 
the dates and time periods of observations for the different instruments, and the altitudes of each 
cloud. 


b. Comparison of the three methods 
A) Illustrative case: flight leg 8 

Figure 8a shows the temporal and altitude variation of the CPL lidar depolarization ratio. It 
indicates that the cloud top altitude varies between 15 and 1 1 km. Around 2217 UTC (red arrow), the 
cloud top altitude is constant, and the depolarization ratio ranges between 30% and 50% within the 
cloud, indicating that it is an ice cloud. Figure 8b shows the GOES visible image at this time together 
with the ER-2 plane flight track. The rise in the cloud-top around 2211 UTC is evident as a shadow 
line in the visible image. At 2217 UTC, the ER-2 is over the ocean, making this flight leg a good 
case study because the surface albedo is low. Figure 8b shows that the plane crosses the edge of a 
cumulonimbus anvil at the end of the flight leg (blue and magenta line segments). For this flight leg, 
the LaRC method classified most of the pixels as ice except for a few pixels that were defined as 
liquid water or clear. The MAS and CALIPSO-like techniques classified all of the pixels as ice 
clouds. This difference is likely due to the large discrepancies in spatial resolutions between the 
LaRC method using GOES data (8-km resolution across and 4-km resolution along the flight leg) 
and CALIPSO-like and MAS methods using MAS data (150 m resolution along the flight leg). The 
patterns in optical thickness seen in the MAS (Fig. 6c) and GOES (Fig. 6d) results are generally 
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consistent over the relatively short length of the flight leg. The ice cloud optical thickness is often 
larger than 3 in anvils; hence for the CALIPSO-like method, the cloud optical thickness is 
determined in the IR by comparing the radiative transfer computations with the measured radiance. 
This was not the case for the SIRTA study where the optical thickness was detennined from lidar 
because clouds were semi-transparent. 

Table 6 summarizes the results obtained with the 3 methods for flight leg 8. Both passive 
remote sensing techniques underestimate the cloud altitude as detected by the lidar: The LaRC 
method locates the cloud between 4.5 km below to 0.5 km above cloud base, whereas the MAS 
method identifies it just below the cloud base. Similarly, T LaRC is -27 K greater than T C AL-iike- Hence, 
both passive methods miss the highest part of the cloud, whereas the CALIPSO-like method misses 
the lower part of the cloud where the lidar cannot penetrate. This result is consistent with the analysis 
of Sherwood et al. (2004) who found that the thennal infrared GOES channel brightness 
temperatures correspond to an altitude 1- 2 km below the physical top of even the densest 
cumulonimbus clouds. 

Generally, there is a good agreement between the CALIPSO-like method and both passive 
methods for the optical thicknesses, although the GOES yields a greater maximum value because its 
large pixels include the bright clouds seen in the lower left-hand corner of the MAS image. These 
clouds were not included in the narrow strip of MAS and CALIPSO-like retrievals. The three 
methods yield mean particle effective radii ranging between 22 and 40 pm with variability between 7 
and 18 pm depending on the method. The CALIPSO-like method yields the smallest mean value of 

r e - 

B) All 14 CRYSTAL-FACE flight legs 

Table 5 summarizes the results for 14 CRYSTAL-FACE flight legs. LaRC results are 
reported only for the pixels identified as ice clouds. For 5 flight legs (all three 23 July flight legs, 
flight leg 7, and flight leg 12), no particle effective radius and optical thickness were retrieved using 
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the CALIPSO-like method because the measured BTD values were outside the range of the model 
calculations. Those 5 flight legs correspond to very large values of x La RC and x M as (Table 5) with a 
maximum value of Tl 3 rc equal to 127 (flight leg 12). For such thick clouds, the BTD’s measured by 
MAS in the CALIPSO-like channels have values ranging between -1 and -3.5 K (depending on the 
spectral interval). Such values cannot be reproduced by simulations, whatever the microphysical 
properties of the crystals. All those cases correspond to large convective systems as seen in the 
GOES visible images (not shown). During four of the five flight legs, the lidar depolarization only 
yields aspect ratio of class I, indicating that the cloud particles are likely to be plates or nearly 
spherical particles. This class was found in seven other flight legs. For those 5 flight legs, Table 5 
shows that the LaRC optical thickness is greater than 1 1 and the MAS one is greater than 7. It means 
that the CALIPSO-like method cannot be applied to clouds thicker than Xl 3 rc = 1 1 or Xmas = 7. This 
limitation is a general limitation of the CALIPSO-like technique; it is not the same as the SIRTA 
cases limitation that was only due to the fact that the SIRTA lidar was at ground whereas the MODIS 
instrument was space-board. 

For two cases (flight legs 13 and 14) that also correspond to very thick clouds, the CALIPSO- 
like method was applied to pairs of channels in order to select a particle model for each BTD, but 
none of the particle models could explain consistently the three different BTD values, hence the 
particle effective radius could not be retrieved. This is probably because the cloud particles in those 
flight legs are not represented in any of the CALIPSO-like particle models. Actually, it could 
correspond to very small non-spherical particles like droxtal that have been observed at the top of 
anvils. Infrared and visible optical properties computation are not yet available for such crystals, but 
will be available for the operational version of this method, i.e. for applying it to the CALIPSO 
instruments. For the 8 cases when r e ,cAL-iike could not be retrieved, however, both passive methods 
produced results. As seen in Figure 9, r e ,MAS greatly exceeds its LaRC counterpart for 5 of those 
cases, while results from the two methods agree within + 1 pm for the other 3 cases. 
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Results from the three different methods are compared in Figs. 10 - 12 for the other seven 
flight legs. Figure 10a plots T C AL-iike as a function of T LaRC . Variability is the difference between the 
cloud top and bottom detected with CPL, whereas the LaRC one is a standard deviation of the cloud 
effective temperature, calculated in the box of pixels used for the retrieval. Figure 10b shows the 
probability density function of the difference between T C AL-iike at the cloud middle height and the 
mean T LaRC . Except for one case, T C AL-iike is always less than T LaRC . Actually, as for flight leg 8 
(previous section), the LaRC method often misses the highest (coldest) part of the cloud (cooler) 
whereas the CALIPSO-like method often misses the lowest (wannest) part of the cloud. Exceptions 
include the two cases when the lidar detected a midlevel cloud below the cirrus cloud. In those two 
cases, the optical thickness of the upper cloud was so small that the LaRC technique only detected 
the lower cloud because T LaRC and the low cloud temperature from CALIPSO-like method are nearly 
identical (not shown). In 45% of cases, the difference between TcAL-iike and T LaRC is smaller than 10 
K. This tendency is confinned in Figure 10c and d that compare the cloud altitudes retrieved from all 
3 methods. In 92% of cases, z C AL-iike is higher than z LaRC , and in one case when the lidar detected two 
layers, the difference is 7 km. The difference between z C AL-iike and z MAS , is smaller: z C AL-iike is larger 
than z M as in only 43% of the cases. Otherwise, z M as is higher than z LaRC all the time, and the 
difference is larger than 6 km in 45% of cases. For the thin clouds, e.g., the three flight legs of July 
26 th 2002, the lidar can detect even the thinnest parts of the cloud, as indicated by the range of 
retrieved optical thicknesses. As seen for the SIRTA cases, the range of optical thicknesses retrieved 
by MAS and LaRC tend to be bounded by a greater minimum value than found by the CALIPSO- 
like approach and the passive method cloud heights are, on the whole, less than the lidar-determined 
values (Figure 10c). 

Figure 1 1 compares the optical thicknesses obtained with the 3 methods. The difference 
between Xmas and TLaRC is less than 1 for only 12% of cases, and Xmas is larger than T La RC by more 
than 2 for 50% of cases. As shown in Table 5, this difference is very large for the flight legs without 
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CALIPSO-like retrievals. Actually, the CALISPO-like optical thickness cannot be retrieved from the 
same technique as for SIRTA cases because CRYSTAL-FACE clouds are very thick, hence the lidar 
does not detect the cloud bottom and it is not possible to calculate any lidar optical thickness. Hence, 
we use the infrared radiances to estimate an infrared optical thickness that is called xIRcAL-iike- It 
cannot be compared directly with the other methods because it is an infrared optical thickness 
whereas Xmas and Xl 3 rc are visible optical thickness. When considering a cloud composed of large 
particles, the ratio between visible and infrared optical thicknesses is about 2. Actually, Figure 11a 
and b show that the plot of xIRcAL-iike as a function of x M as and x La RC is better matched with the x = 
2y line than with x = y line, especially for the smallest optical thickness values. Since the CALIPSO- 
like optical thicknesses for the CRYSTAL-FACE data are based on infrared retrievals, any values 
exceeding 3 or 4 are highly suspect. The MAS retrieval tends to yield greater optical thicknesses than 
the LaRC approach. 

Figure 12 shows the comparisons between the particle effective radii obtained with the three 
methods. In only one of those cases where the 3 retrieval methods can be applied, is the effective 
radius retrieval good (Figure 12a and b). For the other cases, the mean CALIPSO-like retrieval can 
be twice the LaRC and MAS values. The agreement between the CALIPSO-like and LaRC results is 
slightly better than that with the MAS method. On average, for the five cases, r e ,cAL-iike is 20.5 pm, 
compared to 1 1 and 13.5 pm for LaRC and MAS. The differences in r e for the two passive methods 
in Figure 12c (37%) are much like those seen in Figure 9 (33%). Table 5 shows that for those 5 flight 
legs, the lidar depolarization ratios lead to aspect ratios Q in classes I, II, and III, he. Q < 1.1. With 
this knowledge, it is difficult to trust the r e ,LaRC retrieval as it only uses column habits, corresponding 
to Q > 1.1. This could explain the large difference between r e ,cAL-iike and r e ,LaRC and the smaller 
differences between r e ,cAL-iike and r e ,MAS- 

c. Comparison with in-situ measurements 
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The particle effective radius results are compared to In-Situ measurements in Table 7. For the 
sake of sampling, the WB-57 aircraft is considered to be collocated with the ER-2 if the time 
difference between the passages of two aircraft over the coincident point (defined in Sect. 4a) is less 
than 22 minutes. Hence, five flight legs have the WB-57 passing over coincidence points to within 
14-22 minutes of the ER-2 flight path. For example, flight leg #12 is at 19.82 UTC for ER2; the 
coincident point for WB57 is 14 minutes later, i.e. at 20.05 UTC; hence, the averaged value and the 
standard deviation of r^sim (Table 7) is estimated using values between 19.97 and 20.13 UTC, i.e. 
10 minutes around the coincident point. A value of r e j n .si tu is estimated at each time by calculating 
<V> and <A> for the all size distribution, from 0.175 to 800 pm and using Eq. 1. Tests have been 
done (not shown) by considering the size distribution only from 1.5 to 800pm in order to know if 
small particles introduce a bias in the results: it only changes the results by 0.5 pm at the most. 
Furthermore, we explained in Sect. 4a that the habit hypothesis for In-Situ leads to a ± 40% error on 
the r e) in_situ calculation. The value of this uncertainty is also presented in Table 7. Only one of those 5 
flight legs corresponds to a case with a possible CALIPSO-like retrieval. For this flight leg 8, as 
noted earlier, the agreement is good between CALIPSO-like, LaRC, and MAS, and the Table shows 
that it is also in a good agreement with in-situ measurements. Nevertheless, the uncertainty about 
r e , in-situ shows that its value can be from 20 to almost 50 pm, so with a large range. Also, the remote 
observations give a mean effective radius weighted over the vertical extent of the cloud whereas the 
In-Situ observations are collected at a given altitude of the cloud. Moreover, the coincidence time for 
this flight leg is 22 minutes, the greatest for the set and that could lead to collocation problem. For 
other cases (Table 7), the In-Situ particle sizes are still in a very good agreement with the other, in 
particular with MAS method. They still have a larger range because of the uncertainty, and for leg 
12, 13 and 14, they tend to be larger than the LaRC ones. This result is not consistent with the fact 
that In-Situ measurements are often biased by small particles. 

5. Discussion and conclusions 
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This study intends to compare ice cloud properties derived from the so-called CALIPSO-like 
method using collocated passive and active remote-sensing observations with two other well known 
passive imager methods (from the CERES team and the MODIS operational cloud product team). 
The CALIPSO-like method requires lidar measurements and is most suitable for thin ice clouds. 
Nevertheless, this method must be validated in different ice cloud conditions (thin and thick) and at 
different latitudes (mid-latitudes at SIRTA and Tropics in Florida). Figure 13 shows the phase 
retrieval from the 2 passive remote sensing methods when the lidar algorithm detects an ice cloud 
(i.e. for all the selected cases). The “liquid water + ice” occurrence means that liquid water and ice 
are found in the sample of pixels, and not that there is mixed phase. Actually, there is no mixed 
phase category for LaRC and MOD06/MAS retrievals. The agreement between the 3 methods 
concerning the detection of ice is better for tropical ice clouds: 70% of tropical clouds are ice for the 
LaRC method against 50% for mid-latitude clouds, and 84% of tropical clouds are ice for 
MOD06/MAS method against 45% for mid-latitude clouds. For tropical clouds, 30% are interpreted 
as ice and water by the LaRC algorithms and 8% by the MOD06/MAS algorithm. Furthermore, for 
mid-latitude cases, 35% for LaRC algorithms and 30% for MOD06/MAS are classified as clear sky. 
There is no case interpreted as only water by the LaRC algorithms, whereas 5% at mid-latitude and 
8% at tropics are only water for MOD06/MAS. These results confirm the importance of lidar 
observations for thin ice cloud retrievals, most of which are found at mid-latitudes. The CALIPSO- 
like technique is then well suited for this kind of cloud, despite its sensitivity to small changes or 
uncertainties in the brightness temperature of each IR channel used, in particular due to the surface 
emissivity or the calibration (Chiriaco et ah, 2004). 

When the retrieval is possible, the comparison between ice particle effective radii obtained 
with the three methods give encouraging results keeping in mind that the three methods are based on 
different approaches and are using different wavelengths. At mid-latitudes, the lidar is on the ground 
and the Spectroradiometer in space, so the CALIPSO-like method has been applied to clouds with 
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optical thicknesses less than 3, in order to be sure that both instruments are sensing the same cloud. 
Considering such clouds, and also thick enough for being detected by passive remote sensing 
methods, the agreement between the three methods is quite good for r e retrieval, independently of the 
quality of the cloud altitude and temperature retrieval. CALIPSO-like and LaRC methods have the 
best agreement for r e at mid-latitude: in 50% of the cases, the difference between r e cAi -iii<e and r ejLaRC 
is less than 3 pm, and r e>C AL-iike is at least 5pm larger than r ejM OD 06 - In the tropics, CALIPSO-like 
method can only be applied to a few cases, also leading to a relatively good agreement between the 3 
methods: in 40% of cases, the difference between r e! cAL-iike and r e% LaRC or r e ,MAS is less than 3 pm, and 
r e ,CAL-iike seems to overestimate r e LaRC and r e ,MAS in 40% of cases. In the tropics, there are often large 
differences between r e>C AL-iike and the two other methods: those cases correspond to cloud with a 
microphysics that is not taken into account in LaRC and MAS methods (typically plates, as shown by 
the lidar depolarisation ratio). Furthermore, there is a good agreement between r e LaRC and r eiM OD 06 
when the cloud optical thickness is moderate, whereas there is a strong discrepancy in very thick 
clouds for which CALIPSO-like method did not give any results. Hence, it seems that whatever the 
method used, the estimation of the particle effective radius in strongly convective clouds in the 
tropics remains a challenging task. One important conclusion concerning the particle size retrieval is 
the complementary role between the different methods: the lidar method is powerful for thin clouds 
but not for thick clouds, and vice-versa for the passive methods. 

This study contributes to quantify the limits of the CALIPSO-like method. The method can 
contribute to improve the knowledge of thin and subvisible clouds properties where passive methods 
are less reliable. Currently the CALIPSO-like method cannot be applied to derive particle effective 
radius and optical thickness in very thick ice clouds in the tropics because the theoretical radiative 
computation used for this method are unable to reproduce negative BTDs above thick cold ice clouds 
as seen in the observations. Several theoretical scenarios (Chepfer et al. submitted) have been studied 
to explain this discrepancy such as including large variability in the temperature profile around the 
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tropopause level, or different microphysical models, leading to the conclusion that CALIPSO-like 
method using MODIS data fails in presence of nitric acid in ice clouds that is not taken into account 
in the radiative transfer computations. The CALIPSO-like method applied to the real IIR/CALIPSO 
channels that are slightly different of the MODIS ones should give much better results for those 
cases because not influenced by nitric acid absorption. 

Future work will consist in applying the CALIPSO-like method to space-bome observations. This 
will cover a large dataset in various latitude and time periods. The retrieval will be validated against 
in-situ data as well as ground-based retrievals and other satellite retrievals. 
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Table 1: Cloud properties retrieved from 3 methods 
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Table 2: Spectroradiometers and sampling used for each method 



CALIPSO-like method/SIRTA 

LaRC method/SIRTA 

MOD06 method/SIRTA 

instruments 

MODIS/SIRTA lidar 

MODIS 

MODIS 

channels 

8.65, 11.15, 12.05 pm 

0.65, 3.75, 10.8, 12.05 pm 

20 bands from visible to infrared 

sampling 

pixel strip 3-km wide + 1 h around 

pixel strip 3-km wide + 1 h around Terra 

- pixel strip 3-km wide + 1 h around 


Terra overpass 

overpass 

Terra overpass 

- 5x5 pixels area (5x5 km 2 ) centred on 
the SIRTA (for cloud top pressure and 
temperature retrievals) 


CALIPSO-like method/CRYSTAL- 
FACE 

LaRC method/CRYSTAL-FACE 

MAS method/CRYSTAL-FACE 

instruments 

ER-2 MAS/cloud lidar 

GOES-8 imager 

ER-2 MAS 

channels 

8.4, 10.8, 12.0 pm 

0.65,3.9, 10.8, 12 pm 

20 bands from visible to infrared 

sampling 

3x3 pixel array (150X150 nr) nadir 

Mean retrievals from 4x4 4-km pixel 

3x3 pixel array (150X150 m 2 ) nadir 


over length of flight leg 

arrays along ER-2 flight leg 

over length of flight leg 
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Table 3: Mid-latitude ice clouds case. (1): ground temperature from SIRTA instruments, “b ” 
subscript is “ base ”, “T” subscript is "Top”, “e” subscript is “effective”. 


Date 

UTC 

Skin 

Temp 

(K) 

CALIPSO-like 

LaRC 

MOD06 

Comments 

Z b ( km) 
T„ (K) 

Z T (km) 
T t (K) 

^lidar 

z e , (km) 
T e (K) 

tice 

z T , (km) 
T t (K) 

^ice 

5/3/02* 

1110 

281 

7.3 

234 

9.5 

214 

0.8-1 

7.7 

238 

0.8 

9.9 

217 

1.4 

LaRC: 2 ice layers 
MOD06: ice 

2/4/02 

1045 

290 

7 

242 

9.2 

225 

1.5-2 

9.2 

226 

2.9 

10.4 

215 

3 

CALIPSO-like: cloud top not 
detected 

LaRC: 2 ice layers 
MOD06: ice 

1/10/02* 

1055 

290 

8.8 

236 

10.7 

219 

1-4 

8.8 

235 

2.2 

10.1 

225 

2.6 

LaRC: ice and liquid water 
MOD06: ice 

8/10/02* 

1105 

285 

6.1 

251 

11 

214 

1.4-3 

7 

244 

1.8 

9.7 

222 

2.2 

LaRC: ice and liquid water 
MOD06: ice 

14/10/02 

1210 

285 

5.5 

259 

8 

241 

1-3 

10.2 

224 

11 

No cloud 

11 

CALIPSO-like: cloud top not 

detected 

LaRC: ice 

MOD06: ice 

6/11/02* 

1030 

282 

6.9 

247 

10.7 

222 

0.4-1. 2 

6.3 

250 

2.4 

10.9 

217 

3.2 

LaRC: ice and liquid water 
MOD06: ice 

19/12/02* 

1015 

273 

8.6 

230 

12 

202 

0. 8-1.2 

6.8 

239 

0.9 

10.4 

214 

1.3 

LaRC: 2 ice layers 
MOD06: ice 

19/12/02 

1150 

274 

10.3 

216 

11.9 

202 

0. 2-0.3 

6.8 

235 

0.3 

No cloud 

1 

LaRC: 4 ice layers 
MOD06: 

19/2/03 

1025 

271 

8.8 

223 

10.2 

211 

0.2 

No ice 

No ice 

LaRC: liquid water 
MOD06: liquid water 

24/2/03 

1040 

283 

7.5 

234 

8.7 

225 

0-0.1 

2.2 

264 

0.3 

No cloud 

LaRC: ice and liquid water 
MOD06: no cloud 

17/3/03 

1100 

284 

6.8 

238 

7.8 

231 

0.1 

3.6 

257 

0.4 

No cloud 

1 

LaRC: ice and liquid water 
MOD06: ice 

27/3/03 

1140 

291 

9.8 

220 

10.1 

214 

0.06 

No cloud 

No cloud 

LaRC and MODIS: no cloud 

1/4/03* 

1020 

288 

9.4 

224 

10.7 

214 

0.15-0.8 

7.4 

239 

1.3 

No ice 

1.7 

LaRC: ice and liquid water 
MOD06: ice and liquid water 

9/9/03* 

1100 

292 

9.2 

230 

10.1 

224 

0.09-0.1 

6.6 

248 

1.3 

16.4 

? 

1.8 

LaRC: ice and liquid water 
MOD06: ice and liquid water 

15/9/03 

1200 

295 

11.7 

220 

12 

218 

0.03 

No cloud 

No cloud 

LaRC and MODIS: no cloud 

6/11/03 

1140 

287 

10.3 

225 

11.2 

217 

0-0.1 

No cloud 

11.3 

218 

No 

cloud 

LaRC and MODIS: no cloud 

17/11/03* 

1120 

282 

9.3 

227 

11.2 

213 

0.05-0.1 

5.6 

249 

0.6 

No cloud 

2.8 

LaRC: ice and liquid water 
MOD06: ice and liquid water 

9/12/03 

1040 

278 

9.3 

227 

11.2 

213 

0. 1-0.2 

2.5 

261 

0.3 

No cloud 

LaRC: ice and liquid water 
MOD06: no cloud 


* When the comparison of the 3 methods is possible (semi-transparent cloud for CALIPSO-like technique and 
ice cloud found for LaRC and MOD06 techniques) 
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Table 4: Results of the three retrieval methods for 8 October 2002 case. 


method 

altitude (km) 

temperature (K) 

optical thickness 

particle effective 
radius (pm) 

CALIPSO-like 

6.1 - 11 

251-214 

1.4 -3.0 

6-22 

LaRC 

6.9 -8.7 

245 - 231 

1.8 ±0.5 

19 ±4.5 

MOD06 

9.7 (top) 

222 (top) 

2.2 ±0.5 

14-37 
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Table 5: Tropical ice clouds cases during CRYSTAL-FACE. 


Day 

(July 2002) 
Start-stop 
times (TJTC) 

Sfc 
temp 
(K) 1 
Leg # 

CALIPSO-like 

(MAS) 

LaRC (GOES) 

MAS (MAS) 

Comments 

z, 

7 (2) 

z^max 

z b 

(km) 

T, 

r r (2) 
1 max 

T b 

(K) 

tIR <3) 

+shape 

class 

T C (K) 

GOES 

time 

(UTC) 

z, 

z b 

(km) 

T 

z t 

(km) 

T 

23 

18.641 

18.655 

302 

#1 

13.8 

11.9 
11.7 

209 

222 

223 

? 

/ 

224 

18.42 

11.7 

11-12 

6.8 

27.9 

14-35 

12.2 

0-35 

Thick cloud. 
LaRC: ice 
MAS: ice 

23 

18.858 

18.870 

302 

#2 

13.6 

11.7 
11.5 

209 

224 

225 

? 

/ 

225 

18.42 


18.9 

11-22 

12.2 

6-8 

Thick cloud. 
LaRC: ice 
MAS: ice 

23 

19.085 

19.091 

301 

#3 

10.1 

8.79 

8.5 


? 

/ 

250 

18.92 

8.3 

8.3 

4.9 

16.9 

16.9 

10.9 

5-10 

Thick cloud. 
LaRC: ice 
MAS: ice 

26 

18.300 

18.359 

303 

#4 

15.4 

12.5 


0.25-1 

II 

259 

17.9 

7.5 

5-10 

7 

0.87 

0.7-1 

10.9 

0-9 

Thin ice cloud, possible low 
layer 

LaRC: water + ice 
MAS: multi ice 

26 

18.482 

18.500 

303 

#5 

15.4 

12.5 

200 

218 


232 

18.75 

1 1.7 
9-15.7 
10.4 

0.87 
0. 5-1.4 

10.9 

0-17 

Thin ice cloud, possible low 
layer 

LaRC: ice 
MAS: multi-layer 
ice+water+clear 

26 

19.000 

19.079 

303 

#6 

15.4 

12.4 


0.5-2 
I. II. Ill 

254 

18.75 

8.7 
6-13 

7.7 

1.3 

0.5-2 

10.9 

0-8 

Thin ice cloud 
LaRC: ice 

MAS: single- layer ice + 
multi mixed 

28 

21.575 

21.578 

303 

#7 

14.6 

13.4 

13.3 


? 

II 

216 

21.59 

12.8 

12.8 

6.7 

68 

57-71 

14.5 

5-48 

Thick cloud. 
LaRC: ice 
MAS: ice 

28 

22.296 

22.321 

303 

#8 

14.6 

11.6 
11.2 



260 

22.32 

9.8 

8.5 

5-12.7 

1.3 

0.5-1. 9 

10.8 

1-3 

Possible low layer 
LaRC: water/ice/ 
MAS: ice 

28 

22.945 

22.973 

303 

#9 

14.4 

13 

13 

210 

215 

215 

No data 

220 

22.99 

12.4 

11-13 

7.7 

21.5 

7-44 

13.5 

1-46 

Missing MAS data 
LaRC: ice 
MAS: ice 

29 

15.702 

15.710 

304 

#10 

11.9 

10.3 

10.3 

+5.7-6 

222 

233 

233 

+263 

0.5-10 

/ 

260 

15.67 

7.7 

7.2-7 .8 
6.6 

2.9 

2-5 

10.7 

4-45 

CALIPSO-like: + low layer 
LaRC: ice 

MAS: single ice + multi- 
layer mixed 

29 

15.811 

15.823 

304 

#11 

11.9 

10.5 

+5-5.5 

222 

231 

+267 

>2 

/ 

261 

15.75 

7.6 

7.6 

6.3 

3.9 

3.9 

10.7 

2-6 

CALIPSO-like: + low layer 
Ambiguous phase 
LaRC: water 
MAS: water 

29 

19.814 

19.827 

304 

#12 

13.5 
13.1 

12.5 

215 

218 

221 

? 

/ 

219 

19.75 

12.5 

12.5 

6 

128 

128 

14.5 

20-62 

Thick anvil 
LaRC: ice 
MAS: ice 

29 

20.300 

20.322 

304 

#13 

13.5 

12.5 
12.5 

215 

221 

221 

2-10 
I. II. Ill 

250 

20.25 

9.5 
6 -12 
7.8 

2.5 

1-3 

10.7 

3-5 

Thick anvil. No CALIPSO- 
like model match 
LaRC: ice 
MAS: ice 

29 

20.424 

20.440 

304 

#14 

13.5 

12.5 
12.5 

215 

221 

221 

1-10 
I. II. Ill 

248 

20.25 

9.9 

7-11 

8.1 

2.3 

1-4 

10.7 

2-5 

Thick anvil. No CALIPSO- 
like model match 
LaRC: ice 
MAS: ice 


1 temperature of dropsonde lowest layer 

Altitude and temperature of maximum lidar backscattered signal (altitude taken into account in the radiative transfer 
calculation for the CALIPSO-like method) 

3 for CALIPSO-like method, optical thickness is from the radiative transfer calculation compared to the observation 
because clouds are often too thick to retrieve the optical thickness directly from the lidar. 
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Table 6: Results of the three retrieval methods for flight leg 8, mean with range in parentheses. 


, , altitude (km) 

method 

temperature (K) 

optical thickness 

particle effective 
radius (pm) 

CALIPSO-like (11.2-14.6) 

(209.0 - 225.0) 

(0.25-3.0) 

(22- 32) 

LaRC 9.8 (5.3-12.7) 

248.0 (234 -268) 

1.3 (0.5- 1.9) 

27.8 (22-40) 

MAS 10.8 

— 

1-3 

(31-38) 
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Table 7 : Comparison between r e j n _su u averaged over 10 minutes, r e cAL-Uke> r e.LaRC, and r eM Asfor the 5 
flight legs with coincidence between the ER-2 and WB-57 flight paths. 


Flight leg 

coincidence time 
(min) 

l*e, In- Situ 
(gm) 

l*e,CAL-like 

(gm) 

r e,LaRC 

(gm) 

r e,MAS 

(gm) 

7 

+ 20 

25.2 (std 14.2) 
uncertainty: ± 10.1 

— 

32.8-34.1 

22.3 -25.9 

8 

-22 

34.4 (std 32.4) 
uncertainty: ±13.7 

22-32 

19.6-40 

31.1 -38.3 

12 

- 14 

43.1 (std 21.1) 
uncertainty: ± 1 7.2 

— 

31.3-33.8 

23.4 -25.7 

13 

- 14 

35.3 (std 24.2) 
uncertainty: ± 14.1 

— 

14.6-22.6 

32 - 39.9 

14 

-22 

37.5 (std 24.1) 
uncertainty: ±15 

— 

16.4-19.8 

32.6 - 39.7 
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Figure captions 

Figure 1: October 8 th 2002: (a) lidar depolarization ratio, (b) composite image using 0.6-, 1.6-, 11- 
pm channels from LaRC retrieval using MODIS data at 1105 UT, (c) LaRC cloud phase: ‘1’ 
is water cloud, ’2’ is ice cloud, ’3’ is no retrieval, ’4’ is clear, ’5’ is bad data, ’6’ is no 
retrieval water cloud, ’7’ is no retrieval ice cloud. SIRTA area is the center of red square. 

Figure 2: x M odo6 retrieved in the 5x5 pixels area against x M odo6 retrieved in the 3 pixels-large wind 
strip, for SIRTA cases. Errorbars are standard deviations calculated in the sample of pixels. 

Figure 3: SIRTA cloud cases (a) TcAL-like = f(TLaRc), (b) TcAL-like = f(TMOD 06 ), (c) TmOD06 = f(TLaRc), 
(d) probability density function of T CA L-iike- T L3 rc, (e) PDF of T C AL-iik e - T M odo6, (f) PDF of 
Tmodo6- T La RC- Dashed lines are x = y. 

Figure 4: SIRTA cloud cases (a) ZCAL-Iike = f(ZLaRc), (b) ZCAL-Iike = f(ZMOD 06 ), (c) ZMOD 06 = f(ZLaRc), (d) 
probability density function of z C AL-iike - z L3 rc, (e) PDF of z C AL-iike - z M odo6, (0 PDF of 
ZMOD06- ZLaRC- Dashed lines are x = y. 

Figure 5: SIRTA cloud cases: (a) TcAL-like - Tl 3 RC = f(XcAL-like), (b) ZCAL-iike - ZLaRC = f(XcAL-like). 

Figure 6: SIRTA cloud cases (a) x CA L-iike = f(x La Rc), (b) x CA L-iike = f(x M OD 06 ), (c) x LaRC = f(x M oD 06 ), (d) 
probability density function of XcAL-iike - XLaRC, (e) probability density function of XcAL-iike - 
Xmodo6, (f) probability density function of x La RC - x M odo6- Dashed lines are x = y. Errorbars are 
standard deviations. 

Figure 7: SIRTA cloud cases: (a) r e , C AL-iike = f(r e ,LaRc), (b) r e , C AL-iike = f(r e ,MOD06), (c) r e , La RC = 
f(r e ,MOD06), (d) probability density function of r e ,cAL-iike - Xe.LaRC, (e) probability density 
function of r e ,cAL-iike - r e , M OD06, (0 probability density function of r e , La RC - r e , M OD06- Dashed 
lines are x = y. Errorbars are standard deviation. 

Figure 8: July 28 th flight leg from 22:17:47 to 22:19:17. (a) CPL depolarization ratio, (b) GOES 
visible image, (c) optical thickness retrieved from MAS data using MAS method, (d) optical 
thickness derived from GOES-8 using LaRC method. Red arrow is the considered flight leg. 
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Figure 9: CRYSTAL-FACE cases for the ones CALIPSO-like retrieval is not possible: r e , M AS = 
f( r e,LaRc)- Dashed line is x = y. Errorbars are standard deviations. 

Figure 10: CRYSTAL-FACE cloud cases (a) T C AL-iike = iTT l a Rc), (b) probability density function of 
TcAL-like - T LaR c, (c) ZcAL-like = f(ZLaRc) in X, ZcAL-like = f(ZMAs) in O, Z La RC = f(ZMAs) i n V, (d) 
probability density function of z C AL-iike - z LaRC , probability density function of z C AL-iike - z MAS , 
probability density function of z LaRC - z M as- Dashed lines are x = y. 

Figure 11: CRYSTAL-FACE cloud cases (a) xIR CA L-iike = f(x L aRc), (b) xIR CA L-iike = f(x M As), (c) x LaRC 
= f(XMAs), (d) probability density function of xIRcAL-iike - x LaR c, (e) probability density 
function of xIRcAL-iike - Xmas, (f) probability density function of x LaR c - x M as- Dashed lines are 
x = y and dotted lines are x=2y. Errorbars are standard deviation. 

Figure 12: CRYSTAL-FACE cloud cases: (a) r e , C AL-iike = f(r e ,LaRc), (b) r e , C AL-iike = f(r e ,MAs), (c) r e , La RC 
= f(r e ,MAs), (d) probability density function of r e ,cAL-iike - r e , La RC, (e) probability density 
function of r e ,cAL-iike - A, mas, (f) probability density function of r e ,LaRC - r e ,MAS- Dashed lines 
are x = y. Errorbars are standard deviation. 

Figure 13: Thennodynamical phase detected by LaRC (in white) and MOD06/MAS methods when 
CALIPSO-like method leads to ice. 
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Figure 1: October 8 th 2002: (a) lidar depolarization ratio, (b) composite image using 0.6-, 1.6-, 11- 
pm channels from LaRC retrieval using MODIS data at 1105 UT, (c) LaRC cloud phase: ‘1 ’ is water 
cloud, ’2 ’ is ice cloud, ’3 ’ is no retrieval, ’4 ’ is clear, ’5 ’ is bad data, ’6 ’ is no retrieval water cloud, 
’7’ is no retrieval ice cloud. SIRTA area is the center of red square. 
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Figure 2: Tmodo6 retrieved in the 5x5 pixels area against Tmodo6 retrieved in the 3 pixels-large wind 
strip, for SIRTA cases. Errorbars are standard deviations calculated in the sample of pixels. 
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Figure 3 : SIRTA cloud cases (a) T CA L-iike = f(T La Rc), (b) T C AL-uke =f(T M OD06), (c) T MO do6 = f(T La Rc), 
(d) probability density function ofT C AL-hke- T LaRC , (e) PDF of T C AL-uke~ T M odo6, (0 PDF of T MOD 06- 
Tlurc ■ Dashed lines are x = y. 
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Figure 4: SIRTA cloud cases (a) z C AL-uke = f(zLaRc), (b) z C AL-Hke =f(zMOD06), (c) z MOD 06 = /(zlorc), (d) 
probability density function of z CA L-iike- z LaRC , (e) PDF of z C AL-uke- z M odo6, (f) PDF of z MOD 06 — z La RC- 
Dashed lines are x = y. 
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Figure 6: SIRTA cloud cases (a) TcAL-uke =f(HaRc), (b) TcAL-uke =f(T M OD06), (c) T La Rc =f(^MOim), (d) 
probability density function of XcAL-iike - % urc, ( e ) probability density function of TcAL-iike - Fioim, (0 
probability density function of Tlurc - ?modo 6 ■ Dashed lines are x = y. Errorbars are standard 
deviations. 
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Figure 7 : SIRTA cloud cases: (a) r e>C AL-iike = f(r e URc), (b) r c , C AL-uke = f(r eM OD06), (c) r eJLaRC = 
f( r e,MOD 06 . ), (d) probability density function of r e cAL-uke - r e,LaRC> ( e ) probability density function of 
r e , CAL-Hke ~ r e,M 01)06, (f) probability density function of r eLaR c - r eM ooo 6 ■ Dashed lines are x = y. 
Errorbars are standard deviation. 
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Figure 8: July 28 th flight leg from 22:17:47 to 22:19:17. (a) CPL depolarization ratio, (b) GOES 


visible image, (c) optical thickness retrieved from MAS data using MAS method, (d) optical thickness 


derived from GOES-8 using LaRC method. Red arrow is the considered flight leg. 
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Figure 10: CRYSTAL-FACE cloud cases (a) TcAL-iike = f(TiaRc), (b) probability density function of 
T CAL-like ~ T LaRC , (c) Z C AL-like = fClaRc) VI X, Z C AL-!ike = f(?MAs) In O, Z LaRC = f(z M As) In V, (d) probability 
density function of zcAL-uke - ZLaRC, probability density function of zcAL-uke - zmas, probability density 
function of z/ a Rc - zmas ■ Dashed lines are x = y. 
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Figure 11: CRYSTAL-FACE cloud cases (a) TlR CA L-uke =f(F a Rc), (b) TlR CAL -uke = f(^MAs), (c) r La RC = 
f( Fias), (d) probability density function of TlRc A L-nke - Fqrc, (e) probability density function of tIRcal- 
like - Ft as, (0 probability density function of Furc - Fias- Dashed lines are x = y and dotted lines are 
x=2y. Errorbars are standard deviation. 
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Figure 12: CRYSTAL-FACE cloud cases: (a) r e>C AL-uke =f(r e .LaRc), (b) r^cAL-uke =f(r e ,MAs)> ( c ) Clcrc = 
fCe.MAs), (d) probability density function of r e cAL-uke - r e ,LaRC , ( e ) probability density function of r e cAL- 
like - r e,MAS> (f) probability density function of r eLaR c - r e M AS ■ Dashed lines are x = y. Errorbars are 
standard deviation. 
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Figure 13: Thermodynamical phase detected by LaRC (in white) and MOD06/MAS methods when 
CALIPSO-like method leads to ice. 
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